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ABSTRACT 

A general understanding of the effect of longitudinal conduction on 
fluid temperatures for a multistream heat exchanger is obtained from 
a plot of a characteristic function. Mathematically, three corre- 
sponding sets of data are found to be real and distinct except for one 
possible pair of zero-double- roots for the balanced flow case. This 
assures four different forms of solutions. The function suggests the 
importance of the diagonal terms in the exact solutions. Therefore, 
an order of magnitude analysis is made on the basis of this function 
and the diagonal terms. A new descriptive parameter is derived and 
the importance of flowrate can be seen. An inefficiency plot was deter- 
mined from an example analysis of a two-stream symmetric case. 
Performance equations, correction formula, and quantitative design 
charts are presented. The effect of large temperature differences on 
the ends is discussed mathematically. 
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DEFINITION OF SYMBOLS 


Symbol 

Definition 

Units 

A 

Heat transfer area 

ft 2 /ft of length 


Total area subject to longitudinal 
conduction 

ft 2 

a nk 

Defined by equation (55) 1/(1 + Xj^/a n ) 


c 

Specific heat of fluids 

Btu/lb -°F 

C' 0 > C o> 

C^, C - Integration constants 


F(X) 

Characteristic function defined by 
equation (4) 


f(N K ) 

[l/ 2p(l-p)] defined by equation (64) 


g(N K ) 

[1/ 2p(l+p)] defined by equation (65) 


h 

Heat transfer coefficient 

Btu/hr-ft 2 -°F 


^N - l c > Loss of effectiveness due 
to the effect of longitudinal conduction 


K 

Conductivity of wall materials 

Btu/Hr-ft-°F 

k,,k 2 

Conduction heat flux at warm (or 
cold) end 


L 

Length of the exchanger 

ft 

m 

Mass flow rate 

lbm/hr 

n k 

(KA c /mc)tt, Modified longitudinal- 
conductivity-flowrate- ratio, dimensionless 

Nl 

NTU = a(L/2), Modified length p 
parameter, dimensionless 



Ng KA c /(mc) m L, defined by equation (46) 
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DEFINITION OF SYMBOLS (Continued) 


Symbol 

N Mo 

p (n k ) 
Q 

R (N k ) 

s (n k ) 


■n 


• w 


U 

WOLC 


WLC 


Definition 


Units 


(L 2 /2)(hA/KA c ) 

1/2 

(1 + 2/N k ) 


Total heat transferred 
(i-f + g) l/2 
(i-f + g) ,/a 


Btu/hr 


Fluid temperatures, dependent °F 

variables 

Wall temperature, dependent °F 

variable 

Inlet temperature at warm end °F 

Inlet temperature at cold end °F 

Overall heat transfer coefficient Btu/hr- F-ft 2 


The case where the effect of 
longitudinal conduction is neglected 


The case where the effect of 
longitudinal conduction is included 


X 


Independent variable in length 
direction ft 


Y 


Dependent variable in F(\) versus \ 
plot 


Greek Letters 


a n (hA/mc) n , Asymptotic ratios 

p n Roots for the case of WOLC 
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DEFINITION OF SYMBOLS (Concluded) 


Symbol 

^•n 

£ 

e 

Subscripts 

C 

infl 

k 

m 

max 

n 

N 

i 

o , o 
Pi (+) 

qj (-) 


r 


Definition 

roots for the case of WLC 

<T" - Tp) (mc) p /(T'- - T4) (mc) m . 
Effectiveness on the thermodynamics 
base for two-stream exchanger 

(T - T'p)/(T'' - T'q), Dimensionless 
temperature Field 


The case of WLC 
Point of inflection 

Subscripts related to number of roots 
Minimum value 
Maximum value 

Number of streams, or total number of 
streams 

The case of WOLC 
Number of two special roots 

Number of streams flowing in 
positive direction 

Number of streams flowing in negative 
direction 

Matrix rotor in regenerator 


Units 


I 
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A NEW CONCEPT TO THE GENERAL UNDERSTANDING 
OF THE EFFECT OF LONGITUDINAL CONDUCTION 
FOR MULTISTREAM COUNTERFLOW HEAT EXCHANGERS 


By 

C. L. Pan, N. E. Welch and R. R. Head 
George C. Marshall Space Flight Center 
Huntsville, Alabama 


SUMMARY ^ 

3oin 

Considerable interest has recently been given to the design of 
gas flow heat exchangers with large heat transfer areas operating at 
low flow rates. However, these design conditions tend to magnify 
the effect of longitudinal conduction. The mathematical description 
of two- stream and multistream counterflow heat exchangers, including 
longitudinal conduction in the separating walls, by the authors, and 
work of others, have produced exact solutions or approximations that 
are too involved algebraically. Consequently, the relationship between 
the mathematical results and the physical phenomena could not be as 
clearly expressed as with conventional methods. A new concept is 
presented in which a plot of a characteristic function is used to obtain a 
general understanding of the effect of longitudinal conduction for both 
the two-stream and multistream cases. 

The function suggests the importance of the diagonal terms in the 
exact solutions. Therefore, an order of magnitude analysis is made 
on the basis of this function and the diagonal terms. A new descriptive 
parameter was derived and the importance of flowrate can be seen. An 
inefficiency plot was determined from an example analysis of a two- 
stream symmetric case. Performance equations, correction formula, 
and quantitative design charts are presented. The effect of large 
temperature differences on the ends is discussed mathematically. 


INTRODUCTION 


Kays and London [1] pointed out that considerable interest has 
recently been given to the design of gas flow heat exchangers with 
large heat transfer areas operating at low flow rates. Of course, 
these design conditions tend to magnify the effect of longitudinal 
conduction; because of this, the authors have considered the mathe- 
matical description of two- stream (Ref. 2) and multistream (Ref. 3) 
counterflow heat exchangers including longitudinal conduction in the 
separating walls . However, these developments, as well as those 
of previous authors such as Dr. Hahnemann [4}, and Landau and 
Hlinka [5], have produced exact solutions or approximate techniques 
that are too involved algebraically. Consequently, the relationship 
between the mathematical results and the physical phenomena could 
not be expressed as clearly as with the conventional methods. A 
new concept is presented in this report in which a plot of a charac- 
teristic function is utilized to obtain a general understanding of the 
effect of longitudinal conduction for both the two-stream and multi- 
stream cases. A simplified method based on this function is 
presented. 

It is only necessary that (mc) r /(mc) m > 5 for the periodic flow 
to approach the direct transfer type behavior. The analysis and 

physical interpretations presented will aid the reader in the 
analysis and interpretation of periodic flow regenerators for which 
only approximate and numerical results are presently available (Ref. 
6 and 7). 


MATHEMATICAL MODEL 

DIFFERENTIAL EQUATIONS 
THAT DESCRIBE THE SYSTEM 

The differential equations that describe an n- stream heat 
exchanger were derived by Sze and Cimler [8]. They assumed 
infinite conductivity normal to flow for a reversing heat exchanger. 
Application of these equations to plate and fin type heat exchangers, 
including finite conductivity in the direction of flow, is discussed 
in Reference 3. The mathematical description of the two-stream 
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exchanger is, of course, a special case of the general system. The 
differential equations for the multistream case, including longitudinal 
conduction in the separating walls, are written as follows: 


dT 


(j^mc) n + (hA) n (T n - T w ) 0 


n 

S ( i mc) n 
n = 1 


dT 


— = KA C 
dX c 


d 2 T 


w 


dX' 


( 1 ) 


(2) 


The boundary conditions for constant heat flux at the ends of the 
separating walls and constant fluid inlet temperatures are written as 
follows: 


X = 0 

T . - T n 

A pi “ x p 

dT W(o) _ y 

dX 

X = L 

T • - T* 

dT "‘ L) . K 
dX 


(3) 


CHARACTERISTIC FUNCTION 


The following form of the characteristic equation is found to be 
significant. It is derived in Appendix A. 

n 

1 ' (4) 


F(X) =2 


= KA C \ 


n=l X/(hA) +l/(±mc) 
n n 


K = 0 


(5) 


The roots of equation (4) are determined graphically. An example 
is given for a four-stream unbalanced flow heat exchanger. Figure 1 
is the F(X) versus X plot for this case. There are (n + 1) = 5 discon- 
tinuous curves. The intersections of the curves with the X-axis are 
the roots for the WOLC case, and the intersections of the curves with 
Y = KA C X are the roots for the WLC case. The total conductivity 
(KA C ) is the slope of the straight line that shifts the roots from Y = 

0 to Y = KA C X. 

The form of the solution of the set of differential equations 
depends on the type of roots that are obtained from equation (4). 
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CHARACTERISTIC FUNCTION F(\) VERSUS X 
DR FOUR- ST REAM HEAT EXCHANGERS 






T 


1. If all the roots are real, 
n 

= ^ C k exp (\ k X) 
k = 1 


the solution is written as follows: 


( 6 ) 


2. Since the coefficients of the differential equations are all 
real, complex roots appear as conjugate pairs. If only one pair of the 
roots are complex numbers, the form of the solution will be 

T = (C 1 cos \ T X + C 2 sin X^X) exp (X 2 X) + N exp (X 2 X) (7) 

k = 3 


3. If all the roots are real, but one pair are double roots, 

that is 

X,= X 2 = \ 


the solution is 

T = C, exp (XX) + C 2 X exp 


k 

(KX) C k exp (\ k X) 

k = 3 


( 8 ) 


If three roots are equal, the solution is 

T = C, exp (XX) + C 2 X exp (X X) + C 3 X 2 exp (XX) 
k 

+ ^ c k exp (x k X) 
k = 4 


(9) 


4. Any other combination of the following terms 
exp (XX), X exp (XX), X 2 exp (XX), ... cos XX, 
sin \X . . . may appear in the solutions. 


( 10 ) 


Although other mathematical proof can be presented, it is easy 
to explain the possible types of roots from the plot of the characteris- 
tic function. If n values of a n = (hA/mc) n are real, there are (n+ 1) 
curves for which -co ^ F(\) ^ co, -go £ F(k) ^0,0^ F(\) ^ + co, and 
each curve intersects with Y = 0 and Y ~ KA C \ producing n + 1 real 
roots. Therefore, there will be no complex roots, double roots, or 
triple roots, and no solutions in the form of equations (7), (9), and 
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other combinations. But the form of equation (8) is still possible 
since the central curve may pass through the origin, i. e. , one of the 
roots is zero. If \ Q = 0, then from equation (5), = 0, \' Q =\ Q - 0. 

The characteristic equation may, therefore, have one pair of double 
roots equal to zero, which produces a solution of the form of equation 

(8). The simplified form that results is 

k 

T . c , + c !X+ y C k exp (\ k X) (10) 

k = 1 

MATHEMATICAL SIGNIFICANCE OF 

BALANCED FLOW AND UNBALANCED FLOW 


The term "unbalanced flow" for multistream heat exchangers is 
precisely defined as "the condition in which the sum of the time rates 
of heat capacity in the positive direction is not equal to that in the 
negative direction, " that is 


5 > 

n = 1 

It follows 

i<+ 

n = 1 


n 

mc )n + J mc )n * 0 

n = i + 1 

for "balanced flow" that 
n 

mc )n + J (" mc) n = 0 
n = i + 1 


(ID 


( 12 ) 


The central curve of the F(X.) versus ^ plot passes through the origin 
for balanced flow and is slightly shifted for unbalanced flow. This is 
readily illustrated by substituting \ = 0 into the characteristic function 
to obtain 

i n 

F(o> = y <+ mc) n + J (- mc ) n 
n = 1 i + 1 


Therefore, 

F (o) = 0 for "balanced Flow" 

F (o) ^ 0 for "unbalanced flow". 
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Since one pair of zero-double- roots occurs in balanced flow and all 
roots are distinct for unbalanced flow, the basic forms of the solutions 
for these two conditions are different. It is apparent [2] that two 
independent solutions exist. 

TWO SPECIAL ROOTS 

It can be concluded from this discussion that 

1 . for the case of balanced flow and WOLC : 

P'o = Po = 0 

2. for the case of balanced flow and WLC: 

X ’o = = 0 

3. for the case of unbalanced flow and WOLC: 

P'o = 0 

P , >„ * 0 

For the two-stream case: 


UL ! - 1 = Po 

(mcJjj (mc)j 

4. for the case of unbalanced flow and WLC: 


(13) 


X' 0 = 0 

* 0, |x 0 | is smaller than | (3 0 j 

THREE SETS OF ONE- ONE-ONE 
CORRESPONDENT DATA 

If the two roots considered above are excluded, the other roots 
can be traced along each of the F (X) versus X curves (FIG 1). Three 
sets of data, a n , p n , X n , are in one-one-one correspondence with the 
following physical relations: 


( 14 ) 


1. The asymptotic ratios : o- n = (hA/+wc) n 

n •" 1 1 2| • ■ • fi 

If (hAL) n is used to replace (hA) n , to make Q n dimensionless, o-n is 
just the conventional 2 (NTU) n values. From the F(X) versus X plot 
it is found that -«n(s) are locations at F (X) becomes infinite. 

2. The roots for WLC case : X n n = 1,2, .. .n. The X n (s) are 
the values of X for which F (X) intersects the line Y = KA C X. 

3. The roots for WOLC case : p n n = 1 , 2, . . . n. The p n (s) 
are the values of X for which F (X) intersects the Xaxis. It is inter- 
esting to note that for this case: 

P, = -® p n = +® (where n is the last stream) 

If the flowrate and surface conductance of each stream are given, 
these three sets of data can be calculated from the characteristic 
function. 

THE PHYSICAL IMPORTANCE 
OF EXTREME CONDITIONS 

1. When the total conductivity in the direction of flow is 
infinite, longitudinal conduction dominates. 

X n — ■(— <* n ) n = 1, 2, . . .n 

(1 + X n /or n ) 0 (15) 

X Q — 0 for unbalanced flow 

2. When the total conductivity in the direction of flow is finite, 
both conduction and convection are important. 

3. When the conductivity in the direction of flow is zero, 
convection heat transfer dominates. For the WOLC case: 

X n — P n n = 1,2, .. .n 

X, -(-«)= p, (16) 

X n — (+®) = p n (n is the last stream in this term) 

^■o ~ * Po 
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GENERAL SOLUTIONS 


All of the possible roots of the characteristic equation are now 
clear. The general solutions are readily obtained for the following 
cases: 


1. for balanced flow and WLC: 
n 

T = C 1 + C X + V C . exp (\, X) 
w o o wk r k 

k= 1 


T n =C’ 0 + C 0 (X-l/« n )+ V 7 77 7 exp(k k X) 

Cty ! + X. k /or n 


( 17 ) 


(18) 


2. for unbalanced flow and WLC: 

n 

T w = C 'o + C o ex P ( X o X > + 2 C wk ex P <Mc X > * 19 > 

k = 1 


n 


= C 'o + 


1 + 


\o^ a n 


n 

exp (k Q X) + ^ 


""wk 


k= 1 


1 + X.K. / a 


exp (X k X) (20) 


n 


3. for balanced flow and WOLC: 

n- 1 

T w = C 'o + C o X + I C wk ex P (Pk X > 

K = 2 

n " 1 C 

T n ‘ C'o + 0 o ,<X- !/.„>♦ ^ TT “P "V* 

k — 2 

4. for unbalanced flow and WOLC: 

n- 1 

T w = C 'c + C o e *P «>o X > + k I 2 C wk <?k X > 


( 21 ) 


( 22 ) 


(23) 
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( 24 ) 


n- 1 


Tn = c '° + 'Trfir exp tp ° x> * s r:y x > 


A general procedure for the determination of the (n+1) integration 
constants, C* 0 , C Q , C w ^, from the boundary conditions is presented 
later. 


A GENERAL UNDERSTANDING FROM 
A TWO-STREAM SYMMETRIC EXCHANGER 


The temperature fields obtained from the exact solutions are 
necessary for the general understanding of the problem. However, 
the equations are too lengthy for design calculations. Therefore, the 
design techniques are usually based on dimensionless parameters that 
can be used to develop more convenient calculational procedures. 

TWO DESIGN PARAMETERS FOR THE SYMMETRIC 
TWO-STREAM EXCHANGER (h, = h 2 ; m, = m 2 ) 

A dimensionless parameter that includes the effect of longitudinal 
conduction is defined as the modified-longitudinal- conductivity-flowrate 
ratio, which is written as follows: 


KA 

N k = a = 


me 



( 25 ) 


The NTU number may be considered as the modified length, and 
written as follows: 


N t = NTU = - - --- = ( - hA — = a— (26) 

■L* me \ me / 2 2 

The temperature fields, the effectiveness, and the temperature 
differences can be expressed by these two parameters. A formal 
derivation with the physical meaning for Npr will be presented later. 
However, their validity can be clearly seen from the expressions 
presented in the following paragraphs. 
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APPLICATION OF THE 

CHARACTERISTIC FUNCTION 

The characteristic function for the two- stream symmetric case 


F(\) = 


1 


J^ + J- 

hA me 


i^ + J. 


= KAA — i 

c \/a + 1 


hA 


me 


(27) 


-r— = N k (X/a) 

X /a - 1 ^ 


1. For the WOLC case 


1 

P /or + 1 


+ 


1 

P/a - 1 


= 0 


According to equations (21) and (22) 


T w ^ C 'o + C o X 


T = C' + C (X + 1/a) n = 1, 2 (28) 

n o o — ' 

After evaluating the integration constants, the same results as 
Jakob [9] , are obtained. 

2. For the WLC case, equation (27) is used to obtain 


2 X/a = (X/a) (\ 2 /a 2 - 1) 

XL = + 2 N l (1 + 2/N K ) l/2 (29) 

Using equation (17) and (18), the temperature fields are found to 
be similar to those of equations (32), (33) and (34) given in reference 
2 with the exception that they are in terms of N^ and Nj^. 
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MATHEMATICAL EXPLANATION FOR THE 

INCREASE OF END TEMPERATURE DIFFERENCE 
FOR THE WLC CASE 

If t±ie mathematical expressions for the mean temperature 
difference and end temperature difference are considered, an 
interesting observation can be made. The derivation of these 
expressions is given in Appendix B. 


1 + N^ tanh (pN. )/(2pN T ) 

A6 mpan = ^ ^ (30) 

mean N + ^ tanh (pN L )/(2p) + N R /2 + 1 


A0 me an 18 defined by the total heat transferred; that is, 

Q = UAL A6 mean (31) 

If it is compared with the end temperature difference A0(O) or 
A6(L), only one term is found to be different: 

1 + N ^/2 

A0(O) = A0(L) = (32) 

N L + N K 1211111 (pN L )/(2p) + N K /2 + 1 


The temperature field plots in References 2 and 3 suggest that 
longitudinal conduction increases the end temperature difference. 

At low flowrate and high conductivity, a sudden change of fluid 
temperature appears at the exchanger ends. This condition can be 
explained with equations (30) and (32); A0 mean is always smaller than 
A0(O), because the value of tanh (pNjJ is smaller than that of (pNj^), 
and 


tanh (pN L )/(pN L ) < 1 


Therefore, the effect of longitudinal conduction completely 
disappears only when 


tanh (pN L ) = (pN L ) 
for N l (1 + 2/Nj^)’/ 2 = 0 


i. e. 


hA 

me 



0 


(33) 
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Physically, this implies that the total heat transfer surface (or 
the length) is zero, or the flowrate is infinite. 

MATHEMATICAL EXPRESSION OF EFFECTIVENESS 
AND THE RELATION BETWEEN THE WLC AND 
WOLC CASES 


The mathematical expression for the effectiveness of the 
symmetric two- stream case: 


£c = 


N l+ N k tanh (pN L )/(2p) 

N L + N K tanh ( P N L )/(2p) + N k /2 + 1 


(34) 


If equation (33) were true for all values of (pN^ ) and it is substituted 
into equation (34), the following equation is obtained: 


£ 


N " 


N. 


N 


+ N. 


K 


L+ N K N l /2 

n l /2 + n k /2 + 1 


n l ntu 

N l + 1 NTU + 1 


(35) 


Equation (34) is reduced to the conventional form (equation (1)). 


BOUNDARY LINE BETWEEN THE 
WLC AND WOLC CASES 


The derivation given above does not mean that the conventional 
approach would not be true if tanh (pNj^) 4 (pNj^). This can be 
explained by use of an example from the boundary layer theory. The 
velocity in the boundary layer approaches the stream velocity 
asymptotically. However, the boundary layer is thin and is defined 
as the distance from the wall where the velocity differs from the free 
stream velocity by 1 per cent. If this criterion is adopted, a boundary 
line between the WLC and WOLC cases can be found. An algebraic 
expression for this line for the symmetric case is 

N L = 30 N K (1/p + l ) ~ 1 (36) 


This equation approximately describes the values of Nl and Njj for 
which the effectiveness calculated from equation (34) will differ from 
that calculated from equation (35) by less than 1 per cent. It pre- 
scribes the region for which the designer may ignore the effect of 
longitudinal conduction [10] . 
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The 6-NTU chart with Nj^ as a parameter is illustrated in Figure 
2. The curve for Nj^ = 0 is taken from Reference 1, or calculated 
from equation (35). These revised curves are useful for gas flow 
heat exchangers with large heat transfer areas operating at low flow 
rates . 

The £-Nl chart with Nj^ as a parameter and the £-Nj^ chart with 
as a parameter are shown in more detail in Figure 3 and Figure 4. 

BOUNDARIES OF THE 

REGION OF INEFFICIENCY 


Division of equation (30) by equation (32) gives the ratio A8 
A0(O) 


A6 

mean 

A8(0)^ 


N k tanh (pN^)/(2pN L ) + 1 

N k /2 + 1 


(37) 


The designer will prefer to keep this ratio as large as possible. 
Unfortunately, the increase of tanh (pNj__) with respect to (pN L ) 
approaches 1 asymptotically for (pNjJ > 2.7, that is 


(pN L ) > 2.7 0. 99 < tanh (pN^) < 1 


and then 

A6 m.an a N K /2pN L + 1 

A6(0). N„/2 + l <3 ' 

K 

Therefore, the deterioration of thermal performance becomes worse 
for 


N l (1 + 2/N k ) ,/2 >2.7 


(39) 


or 


XL > 5.4 

The dotted line in Figure 3 and Figure 4 of the ^-N^ and £-Nj^ charts 
indicates conditions 

N L (1 + 2/N K )’ /2 < 2.7 (40) 
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It is clear that in this region a large increase of will not cause a 
large decrease of the effectiveness. 

One of the boundaries of the region of inefficiency is defined as 
follows: 

N L (1 + 2/N K ) l/2 4 2.7 (41) 


In this case an abrupt drop of efficiency may occur for high Nj^ 
numbers. 


SIMPLIFIED PERFORMANCE EQUATIONS 


The region specified by equation (39) is not desirable. However, 
it inevitably occurs in practical problems. Fortunately, the mathe- 
matical expressions can be simplified, and the error introduced is 
very small: 


A0 


mean 


A0(O) = 


N K /(2pN L )+ 1 

N K /(2p) + N k /2 + N l + 1 

V 2 + 1 

N R /(2p) + N k / 2+ N l + 1 


N K /(2p) + N l 
N K /( 2p) + N k / 2+ N l + 1 


(42) 

(43) 


(44) 


Comparison of numerical values calculated from the exact 
solution of equation (34) and simplified equation (44) is as follows: 


n k 

n l 

£ exact 

£ simplified 

1.0 

2.7 

.666 

.666 

1.0 

3.0 

.687 

.687 

1.0 

4. 0 

.741 

.741 

1.0 

5.0 

.779 

.779 

10.0 

2.7 

.547 

.548 

10.0 

3.0 

.557 

.558 

10.0 

4.0 

.588 

. 588 

10.0 

5.0 

.6145 

. 6145 
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Further simplification can be obtained for large conduction 
effects: 


1 . 


If N k > 200 


£ = 


n k /2 + N l _ 1 + N b 


n k + n l 


1 + 2N 


B 


Where Ng - (KA c )/(mc) m L Nj^ - 2 Ng Ng 


2 . If and » 1 


*>mean = 1/2N L 


A6(0)' = 1/2 

i = 1/2 


ORDER OF MAGNITUDE ANALYSIS 


( 45 ) 


(46) 


(47) 


It has been shown with the use of equation (39), that XL = 5.4 is 
a special value for the two-stream symmetric case. It is easy to 
generalize this method to unsymmetric, unbalanced flow and multi- 
stream cases. 

First, it can be predicted from the F (X) versus X curve that 
|P n l > l X nl > \ a n\ 

for each of these corresponding values. If a? m is the smallest of these 
asymptotic ratios, and 

<* m L >5.4 or ( N L^m >2 '*^ (48) 

an order of magnitude analysis will lead to simpler expressions. It is 
realized that for the exponential function, 

exp (-X m L) > 200 exp (+X m L) < 0.005 
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The algebraic equations that result will facilitate analytical work for 
large and more complicated heat exchangers. Numerical results are 
also sufficiently accurate if equation (48) is satisfied. 

FOUR-STREAM EXAMPLE 

A four- stream balanced-flow exchanger is considered as an 
example. The two special roots and the one-one-one correspondent 
data for the WLC and WOLC cases are: 

Asymptotic Ratios (<* n ) Roots (\ n ) Roots (p n ) 



Applying the boundary conditions represented by equation (3) to 
equations (17) and (18) for the WLC case, and if equation (48) is 
satisfied, the (n + 1) simultaneous equation for the six constants: 
C' 0 * C Q , C 1 -C 4 are simplified as follows: 



Equations 49, 50, 51, 52, 53, 54 - above. 

Where a^ = 1/(1 + \ k /a n ) (55) 
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Solving C Q from equations (49), (50) and (51) in terms of C ' 0 and 
simultaneously from equations (52), (53) and (54), one readily obtains 
C' 0 . C 0 , C„ C 2 , C 3 , C 4 . 

THE IMPORTANCE OF THE ROLE OF THE 
DIAGONAL TERMS - a kk 

It is concluded from the F(X) versus X plot and equation (15), that 
KA C - ®, (1 + X k /a k ) — 0 (15) 

a kk = 1/(1 + ^k/^k) 03 (56) 


which are the diagonal coefficients in equations (50) through (53). 

This set of coefficients will increase without limit, and when considering 

the characteristic function. 


(mcyimc^ (mc) 2 /(mc) k l 

1 4 ^ k / °1 1 + X k / o 2 1 + 



It is found that when KA C — co, only one term, akk» i n the left hand 
member of the equation remains significant. The other terms in the 
left hand member can be neglected when a kk >100. Equations (49) 
through (54) are further simplified as follows: 


’o 


^2 

0 

0 

1 


r fl 
U O 


' K r" 

1 

a l 1 

0 

0 

0 

-1 /a, 


c, 


T M 

1 p 

1 

0 

a 22 

0 

0 

-l/or 2 


c 2 


T 1 1 

1 p 

1 

0 

0 

a 3 3 eXalj 

0 

(L-1/u 3 ) 


C 3 


T 9 

1 

0 

0 

0 

a 4 4 e>UL 

(L- I/ 0 - 4 ) 


C 4 


T ' q 

0 

0 

0 

X 3 e X3L 

X 4 e K < L 

1 


.C 0 _ 


_ K 2_ 


Equations (57), (58), (59), (60), (61), and (62) - above. 
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Numerical calculations show that the diagonal terms are the 
predominant contributors to the abrupt change of the fluid tempera- 
tures at the inlets. 

CORRECTION FORMULA FOR THE 
TWO-STREAM SYMMETRIC CASE 

If the method used to obtain equations (49) through (54) is used 
for the two- stream symmetric case, the performance equations (42) 
through (44) will be obtained. This approach yields a somewhat 
different physical interpretation that may further simplify the problem. 
Using 

= a 2 ~~ ~ o, Oj — —O', Xj — — X, X 3 — X 

one obtains 


N L - f(N K ) - g(N K ) 
N L - 2 f(N K ) + 1 


where 


(2X/o) (1 - \/a) 


1/ (2p) ( 1 — P ) 


(63) 


(64) 


g = l/(2p) (1 + p) (65) 

It is noted that f is derived from aj^ = 1/(1 + X^/ar^), but g is 
derived from a^ = 1/(1 + X k /ar n ) for n t k. Therefore, f is usually 
more important, and g can be neglected for Nj^ >50. 

A correction formula can be based on the conventional approach 
of increasing the exchanger length to improve the effectiveness. If, 
for constant N^, the loss of effectiveness due to longitudinal conduc- 
tion could be compensated by an increase of Nj^, the same effective- 
ness could be obtained for different lengths. Utilizing equations (35) 
and (63), the following correction formula can be determined. 

N L , C = f + « + N L' N J 1 - f + «> (66) 


for 


N L, N (1 + 2/ N K ) ,/2 > 2 * 7 
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>2.7, equation (66) can be 


If N K > 50, and N L N (1 + 2/N K ) l/2 
further simplified to 

N L, C = f + N L,N < 67 > 

AN EXAMPLE CALCULATION 

USING THE CORRECTION FORMULA 

1. It is of interest to consider the per cent increase of the 
length of a symmetric exchanger that is necessary to maintain the 
same effectiveness calculated from the equation for a WOLC case 

(me) = 8 Btu/hr°F, hA = 160 Btu/hr ft°F, 

KA C = 1 ♦ 0(Btu/hr°F) .ft L = 1 ft 

N l n = (hA/mc) (L/2) = 10 > 2.7 KA c = 1.0, 

N k = (hA/mc) (KA C me) = 2.5 

p = (1 + 2/N k )' /2 = 1.34164 f = 1/(1 - p) (2p) = -1.091 

g = l/2p(l + p) = 0. 159 N l c = f + g + N l N (1- f + g) = 21.67 
AL% = 21.67/10 - 1 = 116.7% 

It appears to be impracticable to double the heat exchanger length in 
order to compensate for the loss due to effect of longitudinal 
conduction. 

2. If the flowrate is doubled in this example, that is,^ 

(me) =16.0 then N, = 5. 0 > 2. 7 

L j IN 

N K = 0.625 p = 2.04939 f= -0.232 g=0.08 N^ c = 6.408 

AL% =28. 16% 

Compensation by increasing the length appears to be more reasonable. 
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3. 


that is (me) = 4 


If the flowrate is reduced to one-half. 


N L n = 20>2 - 7 n k = 10 P= 1.09545 f= -4.7819 
g = 0.218 c = 116.434 AL% = 482% 


Compensation by increasing the length becomes more reasonable. 
Therefore, the length correction method is practical only for the 
higher flowrate range. The 10 per cent correction for this effect as 
mentioned in Reference 8 would be reasonable only for some flow- 
rate higher than 22. 


THE IMPORTANT ROLE OF FLOWRATE AND 
THE DERIVATION OF THE CONDUCTIVITY- 
FLOW RATE-RATIO N R 

The above calculations illustrate that the effect of longitudinal 
conduction is very sensitive to flowrate. There is a mathematical 
derivation that supports this point. Recalling the diagonal terms, 
a^, which are over-weighted due to the effect of longitudinal 
conduction, 

a kk = 1/(1 + X-fc/ok) 

For the two- stream symmetric case 


1 + V"* = 1 - (1 + 2/ V' /! 


By use of the Binomial Theorem, 

— = 1 - |l + 1/2 + LllhlW 


a kk 


N 


K 


(n K ) + 


If N K > 2 


1 

a kk 


= 1 - (1 + i/n k ) = - 


n k 
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Therefore, 


a kk 


-n k = 


(hA) (KA C ) 
(me) 2 



( 68 ) 


Equation (68) illustrates the conclusions that have already been 
observed in this and previous papers. 

1. Flowrate has more effect than any other factor 

2. The new parameter Nj^, which stands for the modified con- 
ductivity-flowrate-ratio, will give a much better grasp of the effect 
of longitudinal conduction. 

3. The parameter Nj^ is proportional to the total conductivity 
and the heat transfer coefficient and inversely proportional to the 
square of the flowrate. 

LOWER BOUNDARY OF THE INEFFICIENCY REGION 

It has been pointed out that the efficiency can be increased by 
increasing N^. Bur, on the contrary, the discussion with equations 
(39), (40), and (41) indicates that if Nj^ is decreased below its 
critical value, the effect of longitudinal conduction will be decreased. 
Using equations (35) ancl (36), the loss of the effectiveness due to the 
effect of longitudinal conduction, I e , is found. 

N t N t - f - g 1-f+g 1 

T - t _ t = tL _ (69) 

6 5 N 5 c 1 + N l l + N L -2f l + N L -2f 1 +N l 

I 6 - is plotted with respect to N^, with Nj^ as a parameter (FIG 5 and 6). 

The loss of effectiveness is serious only in the middle range of N^ for 

constant Nj^. Each Nj^ curve declines at both ends of N^ and reaches 
a maximum at a definite value of N-^. 


LINES OF MINIMUM EFFICIENCY AND THE 

UPPER BOUNDARY OF THE INEFFICIENCY REGION 


Since f and g are functions of Nj^ only, differentiation of I e with 
respect to N^, yields the line of minimum efficiency. 

dl e _ -d + g - f) 1 _ 0 

dN L (1 + N L - 2f) 2 (1 + N L ) 2 


N 


L, 


maximum 


= (1 - 2f - R)/(R - 1) 


(70) 
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Differentiating again 


V = °» s(n k^ = d - f + g) ,/3 

dN L z rs - 

The point of inflection is obtained: 

N L, infl * (1 - 2£ - S)/(S - »> < 71 > 

Of course, at points of inflection, the curves change their direc- 
tions of curvature. The rate of increase of the effectiveness due to 
the increase of will gradually decline beyond these points. This 
line is considered as the upper boundary of the inefficiency region. 

At the same time, equation (41) may be used to define the lower region. 
The boundaries are indicated in Figure 5 and Figure 6. Figure 7 is a 
top view of the region which starts at Nj^ = 0.1. 

A region affected by longitudinal conduction is also indicated in 
Figure 7. Note that equation (36) represents approximately the 
upper limit of the region at which the effectiveness differs within 
one per cent from the ideal WOLC case. As N L (1 + 2/Nj^) 1 / 2 becomes 
smaller, 

tanh [N L (1 + 2/N k )’^]-~ N l (1 + 2/N K ) l/z 

then the effect of longitudinal conduction becomes small and can be 
neglected for 

N L < 0.4 (72) 

This value of can be used as the lower limit of the region for 
larger values of Nj^. 

A COMPARISON OF DATA WITH 
PREVIOUS WORK 

It should be emphasized again that the inefficiency region is 
defined only on the basis of the effect of longitudinal conduction. It 
does not apply to the efficiency of the whole system. As shown in 
Figures 2 through 12 of Reference 1, the effectiveness decreases as 
N^ decreases in the WOLC case. If the effectiveness of the system 
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FIGURE 7 THE REGION AFFECTED BY LONGITUDINAL 
CONDUCTION AND THE INEFFICIENCY REGION 






is plotted with respect to with N| y/ j 0 as a parameter. References 6 
and 4, a maximum is reached. The effectiveness is derived from 
equation (44) as follows: 


N, 


+ + N M.' N L I' 7 ’! 


e = 


N L + V^Mo * 1 + VVn + V /N Mo + 1 


(73) 


for (N l ! + N Mo )' /j >2.7 


The expressions for maximum effectiveness and the points of 
inflection are lengthy. 

The numerical results of the author's solution are compared with 
previous methods in Figure 8. 


RESULTS AND CONCLUSIONS 


1 . The roots of the characteristic equation were graphically 
located and classified into two kinds - the two special roots and the 
three sets of one- one-one corresponding data. 

2. The asymptotic relations and the effect on longitudinal 
conduction were related to the slope of Y = KA C \ 

6 & 8 — — + oo 

r i K n — 

KM a nl 

3. All the roots were found to be real and distinct, except for 
the central curve which passes through the origin in the balanced- 
flow case. 

4. The set of diagonal coefficients, a^* were found to be 
significant when longitudinal conduction was appreciable. The 
diagonal coefficients mathematically produce the sudden change of 
fluid temperature at the inlets. 

The flowrate is found to have more effect than other factors. 


KA — 0, 
c 

KA C ~ «. 
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FIGURE 8 COMPARISON OF DATA WITH REFERENCE 
PAPERS FOR TWO-STREAM SYMMETRIC EXCHANGERS 




2 

5. A new descriptive parameter KA c hA/ (me) was surmised 
mathematically. 

6. A critical value N L = 2.7 was determined. 

7. A simplified method for the general case was determined 
from the order of magnitude analysis. 

For the two- stream symmetric case: 

8. Closed-form algebraic expressions were determined for 
the effectiveness and the mean and end temperature differences. 

9. Length correction formulas were devised. 

10. An inefficiency region was determined mathematically. 
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